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Abstract 

Using the effective theory pNRQCD we determine the potential energy of a color singlet quark-antiquark pair with 
(fixed) distance r in three space-time dimensions at weak coupling (or « 1). The precision of our result reaches 
0(a 3 r 2 ), i.e. NNLO in the multipole expansion, and NNLL in a a/AV expansion, where AV ~ aln(ar). We even 
include all logarithmic terms up to N 4 LL order and compare the outcome to existing lattice data. 
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1. Introduction 

The potential energy of a static quark-antiquark pair in 
the color singlet at short distances r is an essential in- 
gredient in the theoretical description of heavy quarko- 
nium. Its large-distance behavior, which is probed e.g. 
in lattice simulations, indicates confinement. The de- 
termination of the static potential in three space-time 
dimensions (3D) represents an important check of the 
methods used for the four-dimensional (4D) calcula- 
tion Jj]-[3[]. The results can also be applied within ther- 
mal QCD. In this paper we determine the 3D static po- 
tential for or <K 1 using the effective field theory "po- 
tential nonrelativstic QCD" (pNRQCD) flU. Unlike 
conventional perturbative QCD, this effective theory 
framework allows for so-called ultrasoft effects, which 
are crucial for consistent results beyond one loop (0- 

01. We discuss the renormalization group structure of 
the 3D static potential and present recent results up to 
0(a 3 r 2 ) and partly N 4 LL precision. We also compare 
these results to existing lattice data. The work presented 
here is based on Ref. ]9[]. 

2. Theoretical preliminaries 

Because in D dimensions the mass dimension of the 
coupling is [g 2 ] = M 4 ~ D , g 2 r 4 ~ D is a dimensionless (ex- 
pansion) parameter and we have (at least) the follow- 
ing physical scales involved in the problem: l/r (soft), 



V 



l/r x g 2 r 4 D (ultrasoft), g* z o (non-perturbative). 
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In order for perturbation theory at the soft scale and the 
pNRQCD multipole expansion to make sense we de- 
mand g 2 r 4 ~ D <K 1, i.e. weak coupling. 
For D = 3 at short distances we find V ~ g 2 \n(rv), 
where v is the pNRQCD renormalization scale. This 
implies that the ultrasoft expansion parameter g 2 IV ~ 
l/ln(rv) <K 1, if we choose v ~ V. We conclude that 
we can use perturbation theory at the ultrasoft scale V. 
Therefore we formally distinguish between the scales 
V and g 2 . Logarithms from the ultrasoft perturbative 
computation will have the form ln( V/ v) and are rendered 
small, if we set v ~ V. Thus it is legitimate to consider 
the ultrasoft regime as perturbative, i.e. the pNRQCD 
loop expansion makes sense (for sufficiently small r). 
In the following we will use the index "B" to explicitly 
denote bare quantities. Parameters without this index 
are understood to be renormalized in the MS scheme. 
Throughout this paper we will use the notation D = 
3 + 2e. In position space the bare singlet potential can 
be schematically written as 



g%Cn(D)r 



-2(n+ 



(1) 



n=0 



The singlet static energy can be considered to be an ob- 
servable for our purposes. It consists of the potential, 
which is a pNRQCD Wilson coefficient, and an ultra- 
soft contributiory, either bare or renormalized: 



E s (r) = V SJI + 8E U S S B = V s + 6E 1 '/ 



(2) 



1 If one has enough precision also non-perturbative effects at the 
scale g 2 should be included. We will address them in Sec. [5] 
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Figure 1 : One-loop contribution to the octet propagator. The dotted 
line represents the A component of the ultrasoft gluon field. 



The soft contribution V s ,b equals the purely perturbative 
bare static potential, which was computed in Ref. ifioll 
up to two loops, i.e. Oicr'r 1 ). It is IR divergent at this 
order. Using pNRQCD in the static limit the ultrasoft 
contribution can be expressed in a compact form at NLO 
in the multipole expansion (but exact to any order in 
the ultrasoft loop expansion) through the chromoelectric 
correlator. It reads (in the Euclidean) 



6E" 



■Vi 



T F 



l (D-l)N c 



cite 



-IAVb 



x (vac\ gB E a B (t)<f>^(t,0)g B E b E (0)\vac) , 



(3) 



where we have defined AV = V - V s . The concrete 
result for the ultrasoft correction in D dimensions up to 
two loops is given in Ref. [9]. It is known at one loop 
since Ref. | IT ] (see also J3-Q) an d was deduced at two 
loops from the results obtained in Refs. lfl3ill4ll . 

g 2 v 2e 

The renormalized coupling constant a - -j— has in- 
teger mass dimension and does not run in three dimen- 
sions as a consequence of the super-renormalizability of 
the theory. 

The bare potentials Vb in position space also have inte- 
ger mass dimensions and following Ref [15] we define 



V B = V + 6V. 



(4) 



6V will generally depend on a and V. In the MS renor- 
malization scheme it takes the form 



6V 



(5) 



from which we can derive the RGE's for the different 
renormalized potentials V in the usual way. 
In pNRQCD at leading order in the multipole expan- 
sion the singlet field of the quark-antiquark system is 
free, i.e. it does not interact with gluons. It is therefore 
renormalization scale independent at 0(r°). Similarly 
the singlet potential is not renormalized at this order: 



SV S = 0(r 2 ) . 



(6) 



For the octet field the situation is different. Even at lead- 
ing order in the multipole expansion it has a residual 



interaction with ultrasoft gluons. The octet potential re- 
ceives an ultraviolet (UV) divergent correction from the 
one-loop self-energy diagram shown in Fig. [T] 

C A 



7(D 



• a + 0(r A ) . 



(7) 



Higher loop diagrams cannot contribute at 0(r°). This is 
because the potentials must appear perturbatively (with 
positive powers) in the Z's. Since a has positive mass 
dimension, the potentials would appear with negative 
powers in higher loop corrections to Eq. 01, which is 
not allowed by renormalizability. By the very same rea- 
son the octet field does not require renormalization at 
0{r°). With a similar argument and Ref. 11411 we find 
V A/B = 1 + 0(a 2 ). 

Solving the corresponding RGE's and matching to the 
soft tree-level result we find from Eqs. (0 and 0: 

AV(v) = -aC A ln(r ve yE ' 2 Vtt) + 0(r 2 ) + 0(e) . (8) 

We now focus on the renormalization of V s beyond 
0(r°). The singlet potential is IR safe up to soft one- 
loop order. At two soft loops in dimensional regular- 
ization IR poles up to 0(l/e 3 ) appear [10]. The ultra- 
soft computation in pNRQCD, i.e. Eq. 01, yields the 
counterterms 



Z% = r 2 AVa 2 Z (2 ; l) + r 2 a 3 Z% 2) 



Z$j = r 2 AV 2 aZ^ + r 2 AVa 2 Z^ 2) + rV Z^f, (9) 

(10) 
(11) 

Z^=0 V n>3, (12) 

where the explicit expressions for the Z y ' are given 
in Ref. Jgt]. These are the complete 0(r 2 ) results. 
Eq. (TTZl) reflects the super-renormalizability of the the- 



ory. Eqs. (l9lfTTb are obtained as follows: 

Zy ' comes from the one-loop 1/e divergence in Eq. 01 

and fixes also Z {2 ' X) and Zy by RG consistency. Z^ 1 ' 2) 
follows from the remaining 1/e divergence in the ul- 
trasoft two-loop computation for 5V S , once all subdi- 
vergences have been subtracted. From this result we 
derive Z v ' again by RG arguments. Because the re- 
spective divergent parts of the bare quantities in Eq. (0 
have to cancel, we can now also determine the missing 

(1 3) 

term Z;, ' without actually performing the correspond- 



ing ultrasoft three-loop calculation [9]. This is possible 
since Zy, must not contain terms with negative powers 
of AV ~ a ln(rv), which cannot be absorbed by the po- 
tential, cf. Eq. 0. The fact that the resulting Zy' 3) is 
indeed independent of AV is a non-trivial crosscheck of 
both, the soft and the ultrasoft calculations. Thus we 
have found the complete RG structure of V s at 0(r 2 ). 
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3. Results 

From the counterterms determined in the previous sub- 
section, we can derive the complete running of the sin- 
glet static potential at 0(r 2 ). By solving the RG equa- 
tions we obtain 



1 



V s (v) = V s (r;v=-) + Vf G (r;v), 



(13) 



where Vf G (r; v) is the running and V s (r; v = -) is the 
initial matching condition, which we have determined 
using the bare soft data of Ref. [ 10]. The full MS results 
are given in Ref. Q9|] . Note that Eq. (fT3l is the complete 
RG improved static potential (i.e. the soft contribution 
to the static energy) up to 0(r 2 ). Adding the finite parts 
of 6E US at one and two loops after minimal subtraction 
to Eq. dT3l and setting v = A V to resum potentially large 
ultrasoft logarithms we obtain 



E s (r) 



C F a\n(r z v z s jre 7E ) + -C F {1C A 



4rifTf)a r 



+ C F ff J r 2 { T C 2 A ln i (rAV) 



1 r 

4 L A' 



+ 7 Ci(2r £ -l-21n2)ln 2 (rAV) 
\9n 2 - 2 



n f T F (c A (2+^-) + C F (5--)) - {n f T F f n - 



+ C 



48 ' v 2 

2 , 1 ,„2 



2 /13;r 2 1 , 1 . 1 



,( + -yi + -h\ z 2- -y £ (l+41n2) 

Al 384 2 /E 2 4 r 



-(11+ln^)) ln(MV)l + 0(a\ 2 In ) . 



(14) 



This is the full result for the static energy up to 0(a 2 r 2 ) 
and ultrasoft NNLL order expressed as a double expan- 
sion in ar (multipole) and 1 / ln(rAV) (ultrasoft). 
The omitted 0(a 3 r 2 ) terms do not contain logarithms of 
rAV. Eq. ( TPfl i is invariant under a change of v up to 
0(^-). The dependence of E s (r) at tree level on the 
factorization scale v s is related to the 3D relic of the 
4D pole mass renormalon and would cancel, if we add 
twice the heavy quark mass to Eq. (TBI) . In this work 
it will however be of no importance. Finally we would 
like to note that the condition v = AV(v), produces a 
v independent scale that is nonperturbative in a and re- 
sums a certain class of logarithms, see Ref. |0]. 
Since for v = AV(v) the ultrasoft logarithms vanish, we 
only have to add the (RG-scheme dependent) constant 
term 

V,(r;v=-)l v +C F C 2 A r 2 a i c 2 fi (15) 

to Eq. dT4~T > to reach N 3 LL order. C2,o can be computed 
perturbatively, but requires a three-loop pNRQCD com- 
putation which has not been performed yet. 



At even higher orders in the a/AV expansion, non- 
perturbative effects start to contribute. In order to study 
these effects related to loop momenta k ~ a, we inte- 
grate out the AV scale. This means integrating out the 
octet field and ultrasoft gluons. The degrees of free- 
dom left are the singlet field and nonperturbative gluons 
with energy and momentum of order a. The resulting 
Lagrangian, including the leading order nonperturbative 
effects at 0(r 2 ), reads 



Xn P =TAsHido-V s (r)-6E» s )S 



r 

Ay 



TrfS^E-rfs} 



(16) 



for the case without light fermions (nr = 0) to which 
we restrict ourselves in the following □ The coefficient 
of the nonperturbative operator in Eq. ( fTol ) is C np = 1 
at leading order in the expansion. This result is ob- 
tained by matching to a pNRQCD tree-level diagram, 
where two gluons couple to the singlet field at 0{r 2 ). 
The interaction with nonperturbative gluons produces a 
shift of the energy which is proportional to the 3D gluon 
condensate: 



SE"" = - - 

S - B AV B N C (D-1)D 



(17) 



The leading ultraviolet divergence of the gluon conden- 
sate has been calculated in perturbation theory at four 
loops 1 16]. The determination of the finite piece re- 



quires lattice simulations 11171 1 1 811 and a computation 



to change from the lattice to dimensional regulariza- 
tion 11911 . Taking the result (in the Euclidean) from the 
last reference and renormalizing the bare expression in 
Eq. ( TTTb in the MS scheme yields 



8E?(v) = - 



C 3 A C F r 2 a 4 



AV 



384 



v 

C A a i 



-I(ln(167r)+ rB )-i)+2B G 



(18) 



where B' 



.(.sum 



-0.2 ± 0.4(MC) ± 0.4(NSPT). This 



result is of the same order as the ultrasoft four-loop 
contribution, i.e. 0(r 2 a 4 / AV). For v = Ay however, 
the In [ctjj] term in Eq. (fT8l is parametrically dominant 
compared to the latter and we will include it in the nu- 
merical analysis of our results in the next section. 



2 If we were to include light fermions there would also be operators 
oc S'^^S. They could generate corrections to the static energy, due to 
the quark condensate, which are of the same parametric order as the 
purely gluonic ones. 
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4. Comparison to lattice data 
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Figure 2: Plots of the analytic results for the force F(r) = in "ro 

units" with nt = in comparison to SU(2) lattice data (black dots): 
The dotted (LO) and dashed (LO+NLO) curves show the tree-level 
and one-loop results from Ref. [10]. The four other curves include in 
addition the new NNLO order (a 3 r 2 ) contributions in the multipole 
expansion. The labels NNLL and N 3 LL refer to the precision in the 
ultrasoft a/AV expansion. We have plotted our results in the MS as 
well as in the MS scheme in order to make the scheme dependence 
visible and set v = AV. Depending on the scheme we have used the 
values c^q and c¥$ given in the text for the N 3 LL plots, which more- 
over incorporate the leading nonperturbative logarithm of Eq. 1181 . 

Now we would like to see how well the short-distance 
3D lattice data can be reproduced by our theoretical 
prediction for the static singlet energy and, on the 
other hand, extract numerical values for C2.0 in Eq. (TT3T > 
from fits to this data. In Fig. [2] we show our 0(cPr 2 ) 
MS results for the static force F(r) = ^MgzMj up t0 
NNLL (from Eq. (TBI) ) and up to N 3 LL including the 
leading nonperturbative logarithm (from Eq. (fl4l + 
Eq. O + first line of Eq. CED) for n f = and N c = 2 
together with the data points from the SU(2) quenched 
lattice simulation of Ref. |2O]0 We use a — where 
ro = 0.5 fm is the Sommer scaleQ To estimate the theo- 
retical uncertainties we also transformed the MS results 
to the MS schem^l and added the corresponding curves 
as well as the previously known curves at LO and NLO 
precision to the plot. From a fit of the NNLL MS curve 
to the data point at the shortest distance, where we 
expect the best convergence of the perturbative series, 
we determined = -0.04. Transforming this to the 
MS scheme gives = 2.64. 



3 There is also SU(3) lattice data available in Ref. (2lll , but it has 
less points at slightly larger distances. 

4 For the determination of a see Refs. I9 H2H1 . 

5 Note, that this also implies AVms — * AV^jg in me logs of Eq. 1141 
etc. giving rise to a residual scheme dependence from higher orders. 



If we compare the LO, NLO and the NNLO curves 
with ultrasoft NNLL precision, we find a convergent 
pattern, in particular in the short distance limit. Unlike 
the multipole expansion the a/AV expansion does not 
converge well, even for the shortest distances that were 
probed on the lattice. Indeed, already at r/ro — 0.22 we 
have a 0.60 (< ■§f^-) for the ultrasoft expansion 
parameter, which is typically accompanied by the color 
factor Ca- The lack of convergence is in particular 
reflected in the huge scheme dependence of the N 3 LL 
results, i.e. the big separation of the respective MS 
and MS curves in Fig. [2] at larger distances. Therefore 
we do not trust in the values for C2,o given above and 
regard them instead only as a rough order of magnitude 
estimate. To improve on these numbers we would need 
lattice data at much smaller distances. For a more 
detailed numerical analysis of our results see Ref. |@]. 
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